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We study theoretically magnetoresistance of graphene with the short-range disorder. The key 
parameter determining magnetotransport properties - the product of the cyclotron frequency and 
transport scattering time, depends in graphene not only on magnetic field H but also on electron 
energy e : LUcTq oc H/e^ . As a result, "quantum" {tOcTq 3> 1) and "classical" (oJcT, <^ 1) regimes may 
coexist in the same sample at fixed H, giving rise to a strong magnetoresistance. We calculate the 
conductivity tensor within the self-consistent Born approximation focusing on the case of relatively 
high temperature, when Shubnikov de Haas oscillations are suppressed by thermal averaging. We 
demonstrate that both at very low and at very high magnetic field the longitudinal resistivity scales 
as a square root of H: [qxx{II) — Qxxiji)]/ Qxxiji) ~ C\fll, where C is temperature-dependent factor, 
different in the low- and strong-field limits. Furthermore, we predict a non monotonic dependence 
of the Hall coefficient both on magnetic field and on the electron concentration. Finally, we discuss 
the case of the charged impurity potential and also find a square-root low-field dependence of 
magnetoresistance near the Dirac point. 

PACS numbers: 72.80.Vp, 75.47.-m, 73.43.Qt 



I. INTRODUCTION 

Study of magnetotransport in low dimensional systems 
is a powerful tool to probe the nature of disorder and ex- 
tract information about localization phenomena. In par- 
ticular, measurements of magnetoresistance of 2D elec- 
tron gas in conventional semiconductors, like GaAs, al- 
lowed one to identify a variety of different regimes, such 
as Drude-Boltzmann quasiclassical transport, weak local- 
ization regime and the Quantum Hall Effect (see Refs. [l| 
and 13 for review). 

One of the simplest theoretical approaches to the 
problem, so called self-consistent Born approximation 
(SCBA), was developed in Rcfs.^--i6i for two-dimensional 
(2D) electrons with quadratic spectrum mostly for the 
case of short-range disorder. SCBA approach ignores 
localization effects. This implies that the relevant en- 
ergy scale (which is the maximum of the chemical poten- 
tial and temperature T) is large compared to the in- 
verse transport scattering time r^^, which coincides with 
the quantum scattering time Tq for short-range disorder. 
The key parameter of SCBA is ujcTq , where Wc is the cy- 
clotron frequency. For weak magnetic fields, i^c^q <C 1, 
calculations^ reproduce classical Drude-Boltzmann re- 
sult. In the opposite limit, ujc^q 3> 1, semiclassical 
Drude-Boltzmann approach fails, and the conductivity 
is given by a sum of contributions coming from well sep- 
arated Landau levels (LL)^"— 

Low-temperature magnetotransport in graphene was 
studied theoretically in Refs. [vl-fl^. In Refs.l3@ a general 
expression for the SCBA conductivity tensor of graphene 
with short-range disorder was derived and analyzed in de- 
tail for the case of well separated Landau levels. Other 
types of disorder were also discussed,—"— including dis- 
order potentials having special types of symmetriesJ^ 
High-temperature magnetoresistance was calculated in 



Ref. [T3I wi thin semiclassical Drude-Boltzman approach. 
In Ref. [ij, a T-dependent interaction-induced contribu- 
tion to the magnetoresistance (on top of a substantial 
positive T-independent magnetoresistance) was observed 
experimentally and analyzed theoretically. 

A specific property of graphene compared to conven- 
tional 2D semiconductors is the linear energy dispersion 
of carriers, 
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resulting in the density of states, which increases away 
from the Dirac point: 
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Here v — 10® cm/c is the Fermi velocity, e is the en- 
ergy counted from the Dirac point and = 2 • 2 — 4 
is the spin-valley degeneracy. Corresponding wave func- 
tions are given by exp(ikr)|xi^), where IXf) is the spinor 
with the components (e""^/^, ±e*'^/2)/%/2, and ip de- 
notes the polar vector of the momentum k. 

Important consequence of the linear dispersion of 
graphene is that the cyclotron frequency turns out to 
be energy-dependent 
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for e > f2. 
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Here H is the magnetic field, m(e) = e/v'^ is the energy- 
dependent cyclotron mass. 
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is proportional to the distance between lowest Landau 
levels [see Eq. ([TT)) below] and I = ^Jeh/Hc is the mag- 
netic length. Another consequence is the energy depen- 
dence of quantum and momentum relaxation times for 
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the short-range scattering potential. Here we by defini- 
tion assume that the random potential is a short-range 
one if its radius R satisfy the following inequalities 



a < i?< A, 
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where a is the lattice constant and A is the electron 
wavelen gth (t his definition is different from one chosen 
in Refs. fl6lfl7l ). Such potential does not mix valleys and 
can be written asiii^ 



V'(r) = uo^(5(r - rj 
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where summation is taken over impurity positions. The 
correlation function of V{r) has the same form as for 2D 
electrons with the quadratic spectrum: {V{r)V{r')) = 
K(5(r — r'). Here k = UiU^ and rii is the impurity con- 
centration. Calculating by golden rule the quantum and 
transport scattering times we find that they are different 
and both energy-dependent i^iii^ 
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Ttr{e) = 2Tq{e) 
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where 7 = It^v^ j k. The difference between tj^ and 
Tq is due to weak anisotropy of the scattering arising 
from spinor nature of the wave functions. Indeed, the 
squared scattering matrix element, |C/k'kP is propor- 
tional to Kxv'IXip)!^ and, therefore, depends on the scat- 
tering angle if ~ ip' . Below we assume that 7 3> 1 and , 
consequently, eTq{e)/h ^ 1. The letter inequality allows 
us to neglect localization effects. We also assume that 
disorder does not affect density of states. This condition 
is also satisfied provided that 7 is large, with an exception 
for exponentially small energies 

bandwidth), which are irrelevant for this paper. Under 
such assumptions the conductivity in zero magnetic field 
is given by the Drude formula: 
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, for H^O. 
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We see that both Wc and Tq are energy-dependent and, 
therefore, the parameter 



^ci q — n 
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can be small or large at the same sample for different e. 
Here we used Eqs. ([3|) and (O and introduced the energy 

= h£L^i > hVL, (10) 

which scales as a square root of the magnetic field: 

e* cx (11) 

As seen from Eq. © , at sufficiently high T the tempera- 
ture window might include both the "quantum" (x > 1) 



and the "classical" (x < 1) regions, so that the total con- 
ductivity might show some peculiarities specific both for 
quantum and classical transport. 

In this paper we calculate magnetoresistance (MR) 
of graphene with short-range disorder assuming that 
max(T, \i) 3> and, consequently, the number of filled 
LLs is large. 

The most interesting finding is related to the case 
e» {H 0). We demonstrate that the dominant 
contribution to the low-field resistivity comes from the 
energy scale e*, which is deep below max(T,/i) in this 
case. This contribution is calculated within SCBA based 
on the approach of Ref. 0. Assuming that temperature 
is not too small 
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(this condition ensures that Shubnikov de Haas oscilla- 
tions are suppressed by energy averaging within the tem- 
perature window) we find that for low fields (e* <C T) the 
relative longitudinal resistivity scales as a square root of 
H: 



0.784 £« 



fe(0) ^ Tcosh2(^/2T) 



H. 
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Here ls.gxx = QxxiH) - fex(O), and Qxx{Q) = TTh/2e'^"f. 
For T jj., MR is exponentially small because the en- 
ergy scale e ~ £* is well beyond the temperature win- 
dow. However, for T > n, low-field MR is quite large 
and increases with decreasing the temperature. From the 
side of the lowest fields, the square root dependence 
is limited by exponentially small fields corresponding to 
- w Ae-^'t/'^. Calculation of MR at < e'J"" 
is not controlled because at such energies 7 is renormal- 
ized to unity,— so that impurity potential becomes ef- 
fectively strong and SCBA fails. One may expect that 
for ^ £™™ MR becomes parabolic which implies that 
one should replace factor e^/T with (e*/r)(e*/e™"')^ m 

Eq. uni). 

We show that the square-root dependence of MR is 
also obtained in the opposite limit of large field, e* — > 00 
(H^oo): 

Aqxx _ f 0.96 e,T/^2^ for T > /i, and > T'^/fi, 
QxxiO) ^ \ 0.68 e*//^, for ^i:$>T, and e, > ^l. 

(14) 

Further, we discuss the behavior of the Hall coefficient 
i? as a function of the magnetic field and the electron con- 
centration and demonstrate that it is a non-monotonic 
function of both variables. 

In Fig. 1 we plotted schematically the dependence of 
the longitudinal resistivity and the Hall coefficient on the 
magnetic field for ^ <^ T. Importantly, square-root MR 
is predicted both for very law and very high fields. More 
detailed pictures are presented below in Figs. 4, 5, and 
6. It turns out more convenient to plot all dependencies 
not as functions of H but as function of e* cx y/H. This 
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Figure 1: Schematic plot of the dependence of the longitudi- 
nal resistivity and the Hall coefficient on the magnetic field for 
fi <^T. Region of exponentially weak fields corresponding to 
£. < Emin is marlted by grey color. Our theory is apphcable 
for higher fields. Vertical dashed lines correspond to e, ~ T 
and e* ~ T"^ / n- 

allows us to present the results for yu <C T and yu ^ T in 
a similar way. 

Finally, we analyze the case of charged impurities and 
show that at the Dirac point the low-field MR also scales 
as y/H provided that static screening of such impurities 
is taken into account. 



value of the numerical coefficient in the low-field asymp- 
totic of MR [see discussion after Eq. ([5^]. Indeed, purely 
classical Drude-Boltzmann approach is valid for e ^ e* 
and fails at relevant energies e ^ where Landau levels 
start to separate. On the other hand, for e <^ e.^ Landau 
levels are well separated and the longitudinal conduc- 
tivity contains the density of states squared. This af- 
ter thermal averaging also leads to a \/^— contribution 
to MR which comes from separated Landau levels (see, 
e.g., Refs.[T8land[2ll) and has essentially quantum nature 
(despite temperature is higher than inter- level distance) . 
In the calculations below we use a rigorous approach 
based on SCBA, which treats both, classical and quan- 
tum, mechanisms of the square-root-magnetoresistance 
on equal footing and allows one to describe crossover be- 
tween classical and quantum regions at e ~ e,. 

B. Self-consistent Born Approximation for 
graphene with short-range disorder 

Inequality ^ ensures that disorder does not mix val- 
leys K and K', so that one may calculate the conductivity 
in one valley and then simply multiply the obtained re- 
sult by the factor 2. The single- valley Hamiltonian of an 
electron in graphene in the perpendicular magnetic field 
reads ji^ii^ 

= 7?o + l^mp(r). 



II. BASIC EQUATIONS 

A. Qualitative analysis 

In the beginning of this section, before turning to the 
rigorous calculations, it is instructive to make some qual- 
itative estimates clarifying the physics of the predicted 
square-root MR. To this end we note that behavior sim- 
ilar to given by Eq. (jlSp may be obtained already within 
classical Drude-Boltzmann. Indeed, this approach yields 
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nh 1 + [wc(e)Ttr-(e)]^ 
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(15) 
(16) 



where we used Eq. Q. Importantly, the second term 
in the r.h.s. of Eq. (|16p is peaked near e = within 
the width on the order of Consider for simplicity 
the case fi ~ assuming that <C T. Averaging of 
Eq. (|16|) over energy within the temperature window 
yields two terms: field independent contribution 2e'^^ /irh 
and the contribution of the peak which is on the order 
of — (2e^7/7rS)(e*/T). Taking into account that at low 
fields transverse conductivity is small and, consequently, 
Qxx « l/o-a;a;, we find AQxx/ Qxx{Q) ^ £*/T (X -/H. It 
turns out, however, that such an analysis yields incorrect 



Hc) = V h 



where e is the absolute value of electron charge, Viinp(i") 
is the random impurity potential. The eigenenergies and 
eigenfunctions of Hq are given by 

e„ = nn sign(n)y2H, ipnA^^v) = e^^^^^ Xn{x - kl^) , 

(17) 

where is given by Eq. ^ and 



Xn{x) = 



h\n\-l{x) 

sign(n) /i|„|(.t) 


ha(x) 



, forn = ±l,±2,.. 
for n = 0. 



(18) 

Here h„{x) are the normalized wave functions of the har- 
monic oscillator with the frequency D, and mass h/sl. As 
seen from Eq. ([T7)) . the energy-dependent cyclotron fre- 
quency is connected with fl by Eq. ^ while the relevant 
energy scale £* is given by Eq. ([T0| and corresponds to 
high LL, so that we can use the SCBA for calculation of 
the resistivity. 

In the SCBA, the electron Green function in the short- 
range potential is given byi 



G(e) = 



1 



(19) 
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where self-energy is found from the following equation 

1 



(20) 



As shown in Ref. 0, S is 2 x 2 matrix having nonzero 
matrix elements between Xn and x_„. However, at high 
energies corresponding to high Landau levels, S becomes 
simply proportional to the unit matrix 

E(e) « I](e) J J ^ , iore:s>nn, 

where 



I](e) = A(£)±ir(e) 



(21) 



(signs + and — corresponds to advanced and retarded 
self- energies, respectively). In this approximation, S 
obeysi 



1 



n=0 ^ 



£ - S 



(22) 



Here Nmax is the ultraviolet cutoff. The sum entering 
Eq. ((22)) can be calculated by using the identity 



1 



^ N-W 

n=0 



In ( ) -7rcot(7rM/'), (23) 



vaUd for Nmax > Re(VK) > 1, Im(VK) < 
Re(l^), Im(W) > 0. From Eqs. ^ and ^ we obtain 
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In 
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TT cot 



S)= 



2(;if7)2 



. (24) 

The logarithm entering Eq. (j24[) is a smooth function of 
e and leads to a linear in e correction to S , which is irrel- 
evant provided that 7 ^ In Nmax (see also discussion of 
SCBA in graphene in Ref. [l^ . We substract this correc- 
tion from S and for simplicity use the same notation S 
for thus redefined self-energy. Then we find from Eq. (p4)) 
the system of coupled equations for A and T : 



To sin[27r(e - A)/fiijJ 



cosh [I-kT /fvjjc] — cos [27r(£ — A)/hLOc 

Tq sinh [27rr/?iajc] 
cosh [ZTTT/huc] — cos [27r(e — A)/hu!c 



(25) 
(26) 



where ro(e) — h/2Tq{e). We also normalize F by its value 
at zero magnetic field introducing the quantity 



z{s) = 



Tie) _ j.(£) 



(27) 



Here i'{e) is the density of states in the magnetic field. 
The solution of Eqs. and ([^ can be found an- 

alytically in the limiting cases e ^ e» (x ^ 1) and 



e < (x > 1) : 



2Fne-"/" sin ( ?^') , for e > 



A 



F„ 



(28) 



-, for e ^ e* 



1 + 2ae-^/^ + 2(2a2 - 1) - — j e"^''/^, for e > e*; 

(29) 



Here 



a = cos 



27re 



(30) 



i9(y) is equal to unity (zero) for \y\ < 1 (|y| > 1) and 



F„ F(e„) = 



/ 2cJc(£„) 



(31) 



In the first line of Eq. we expanded z up to the 
second order with respect to exp(— tt/x). As seen from 
Eqs. ©j© and ([HT]) . F„ actually does not depend on 
n for the case of short-range disorder which we concern 
with: F„ = Fq = const. Here 



(32) 



Fn = — . 



Hence, the widths of different LLs are equal, but the dis- 
tance hujc between neighbor LLs decreases with increas- 
ing e. We also see that self-energy changes periodically 
with e on the scale hujc- Both the period ad the shape of 
these oscillations slowly depends on energy due to energy 
dependence of x. The density of electron states, i^(e), fol- 
lowing from Eqs. (gT]) and (HH) is plotted in Fig. 2 for 
^ ma.x{T, fi). At low energies, a; ^ 1, LLs are well 
separated, while for x S> 1 density of states is given by 
zero-field density, Eq. ([2]), up to exponentially small os- 
cillating terms. 



III. CALCULATION OF THE CONDUCTIVITY 

The conductivity tensor is given by thermal averaging 
of energy-dependent tensor (Ty (e) 



de 



dnp{e) 



cry(e). 



(33) 



The longitudinal conductivity (Jxx{^) is calculated by 
summation of ladder diagramsii The result is given by 
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Figure 2: Density of electron states in graphene with short 
range disorder (solid line) and the derivative of the Fermi 
function for the temperatures T <^ jj, and T ^ ^ (dashed 
lines). 



Eq. (4.13) of Ref. 0. Using Eq. 
result in terms of zle) 



As) 



26^7 



one may rewrite this 

z{ef 



(34) 



For high energies a; — > 0, so that z — )• 1 [see Eq. ([29| ] 
and we obtain Drude result, Eq. (|15|) . For x 3> 1, we find 
from Eqs. p4p and the conductivity near n-th LlJ- 




for le 



<rn, 



(35) 



We note that Eq. (f34| may be obtained from Eq. p5|) 
by replacement of l/Ttr{e) with z{e)/Ttr{s). Hence, the 
only difference of the SCBA result compared to the Drude 
one is the renormalization of the density of states given 
by Eq. (f27l). 

The calculation of axy{s) is more subtle. Simplest 
approximation based on summation of ladder diagrams 
leads to Drude-like formula with renormalized scattering 
time. 



I ( ] - 2e^7 ujc{e)Ttr{e)z{e) 
'^-y^^>^ nh z2(£) + [w,(e)Tt,(e)]2^ 



(36) 



in a full analogy with Eq. In fact, there also exists 

another contribution to transverse conductivity, which 
is expressed via thermodynamical properties of electron 



gas 
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ai'y{e) = ec[dn/dH],. 



(37) 



Here n = n{e,H) is the electron concentration in mag- 
netic field for zero temperature and chemical potential 
coinciding with e. The derivative over H is taken for fixed 
e. Hence, 



(38) 



In fact, crilie) yields essential contribution to (Txyis) for 
e <C e* when Landau levels are well separated (in the 



opposite case of overlapping Landau levels (T^^(e) is ex- 
ponentially small). 

Let us now do the integral in Eq. (p3| . As follows 
from Eq. (|29p . the dependence <Tij{e) contains fast oscil- 
lations on the scale huJc, the shape and the period of the 
oscillations being energy-dependent due to energy depen- 
dence of X. Therefore, for relatively high temperature the 
integration in Eq. p3p may be performed in two steps. 
First, we average aij (e) in Eq. ((33|) by the energy interval 
Se, such that e ^ 6e ^ fujjc{e). Such an averaging "fil- 
ters" the Shubnikov oscillations and results in a smooth 
function aij{e). One can also show that after averag- 
ing a^y can be neglected^^ (provided that 7 1) both 
for overlapping^^ (e e,) and for separated (e <C £*) 
Landau levels. Hence below we use 

axy{e) « aiyie). 

It is convenient to write aij in the following form 

a,j{e) = afJe)ri,j{e) , (39) 



where r]ij{£) are dimensionless factors. From Eqs. ((29 
, (|34p , and ([55)1 we find asymptotical behavior of 77^ 



24e 



Vxx 



ll-24a;2e-2'^/^ = 1 
I r- 

\Cl^/x = Ci£*/£, for £ < e* 



-2-'/^*, for£»£. 



(40) 



1 + 2e-2''/^ = 1 + 2e-2'^"'/^- , for e » e*, 
1 , for e ^ e*. 

(41) 

Here Ci = 8V2/37rV^ ~ 0.68. For arbitrary values of e 
the factors 77.^ have been calculated numerically. The 
dependencies d^xxi^) and ffxyis) are plotted schemati- 
cally in Fig. 3. In this picture we took into account 
that axxis) = <yxx{~£),<yxy{£) = -crxyi~£) because of 
the particle-hole symmetry. In the high-energy asymp- 
totics we keep exponentially small terms [proportional to 
exp(— 27r/a;)], since they are important in the calculation 
of MR in the limit of low temperature and low field (see 
below). It is worth also noting that averaged longitudinal 
conductivity is not exactly zero at e = but saturates at 
£ ~ HQ, at a quite small value: UxxT^^/'^^'^I ~ 1/^7 1- 

As a second step, we calculate dij by replacing (Tij{e) 



with aij (e) in Eq. (|33 



j de 


dnp (e) 







cry(£). 



(42) 



and substitute thus obtained Gij into expression for lon- 
gitudinal resistivity 



Qxx 



(43) 



'xy 



The result of calculations depends on relation between 
tree relevant energies T, ji and e* . The magnetic field 
dependence of Qxx is encoded in the square-root scaling 
of e*. Below we discuss separately the cases of low and 
high temperatures. 
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X >1 




Figure 3: Dependence of g^x and Oxy on energy. The region 
a; > 1 is marked by grey color. 



A. Low temperatures, T <^ fi. 

We start with discussing of high-field limit, ^ /.j. 
Since integral over energy in Eq. (j42p is concentrated in 
the narrow temperature window near e = ^ we can use 
low energy-asymptotics for axx and a^y, see Fig. 3. Doing 
so, and replacing in Eq. (|^^ [—dnp {e)/de] with (5(e — /i), 
we find 



'xy 



26^7 

tt/i 2e2' 



(44) 



We see that (ixx <C cTxy. Therefore, Qxx ~ '^xx/o'xy and 



%/^, for e, > /i. 



(45) 



Next we consider the opposite case £* <C H- In this 
case, there are different competing contributions to mag- 
netoresistance. First contribution is obtained quite anal- 
ogously to high-field limit by replacing derivative from 
Fermi distribution with delta function. In this approx- 
imation, axx = o'xxifJ')-, o'xy = ffxyifJ-)- As follows from 



Eqs. (|40|) and (|4T|) . in the limit of large /i these conduc- 
tivities differ from Drude values by exponentially small 
terms only. It is well known that in the Drude-Boltzmann 
approximation MR is absent. Hence, MR should be expo- 
nentially small. Indeed, substituting (7xxip) and axyifJ-) 
into Eq. (gS]), using Eqs.dlHl), (001), and (gt]), and keeping 
terms on the order of exp(— 27r/x') we obtain 



QxxiO) 



Set 



for < A*- 



(46) 



There arc two other contributions which may compete 
with this exponentially small result. Both contributions 
arise due to finite value of temperature, which was in fact 
assumed to be zero while deriving Eq. (j46p . Let us now 
take into account that the function 



dnp{e) 
de 



1 



4Tcosh^[(e-//)/2r] 



(47) 



is peaked near e = fi within a finite width on the order 
of T. First of all, there exists a correction to magnetore- 
sistance due to a small variation of x{e) within a tem- 
perature window. To find this correction we put rjij sa 1, 
expand aij (e) near e = /j, up to the second order with re- 
spect to e—fi, calculate integral Eq. (j42l) and use Eq. (|43|) . 
As a result, we get a quadratic-in-iJ MR 



^exx 

Pxx{0) 



'2 =.4 



(48) 



where C2 ~ IGtt^/Q. This correction becomes lager than 
Eq. (HS)) for relatively weak fields such that ^ £1, 
where 



(49) 



27r/i 



\u"\s^Ji■'|C2T^)' 



At very low magnetic fields, another contribution comes 
into play, namely, the contribution to integral Eq. (j42p 
from the energies e ^ e* which are well beyond the tem- 
perature window. To find this contribution we first write 
^xx{e) in the following way 



<^xx{e.) 



2e^7 



(50) 



The function f(jj) has maximum at y = and decays 
when y becomes larger than 1, or, equivalently, e becomes 
larger than (see Fig. 3) 



1 - Ciy^/4, for 2/ < 1 (e < £,). 



(51) 



Here we neglected exponentially small corrections to / at 
large e and introduced dimensionless variable y ~ e/e* = 
1/^/x. Let us now assume that s* is smaller than T. 
Then one can replace [—dnF{e)/de] with its value at zero 
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energy ~n'p{0), while calculating the contribution to (Tx 
coming from the region e ^ e*. Doing so, we obtain 



2e^7 



(52) 



where dimensionless constant C3 is given by the integral 



dy f{y) = 1.568. 



(53) 



It is instructive to compare the obtained result with the 
"purely classical" calculation which ignores existence of 
LLs. On the formal level, such an approximation implies 
replacement rjij with unity in Eq. (j42p . Simple calcula- 
tion (see also discussion in Scc lII A[) shows that thus cal- 
culated conductivity can be presented in the same form 
as Eq. ([5^ . where one should replace C3 with 



C?- I dv[l--^\= [dy- ^ 



1 + Ax 



7r/2 = 1.571. 
(54) 



This value is close but different from C3 , so that discrete- 
ness of LLs should be taken into account. Now we take 
into account that axx ^ (^xy, so that Qxx ~ l/cn and 



Pxx{0) 



2C3£.|n^(0)| oc VH 



(55) 



Since T is by assumption much smaller than we con- 
clude that MR is exponentially small: 



^Qxx 
Pxxi^) 



2C3£, 

T 



(56) 



Comparing Eq. ((56)) with Eq. (|48|) we see that the former 
contribution dominates when ^ £2, where 



e2 - ^e-^"^ 



< T. 



(57) 



Therefore, in accordance with our assumption, <C T 
in the regime, when Eq. ([55)1 dominates. 

Looking now more attentively at the above derivation 
one concludes that Eq. ([55]) is valid at arbitrary relation 
between T and ^ provided that £* (but e* > er™)- 
Hence, the low-field MR is given by 



(736* 



^Qxx _ 

Pxx{Q) ^ 2rcosh^(/i/2T) 



oc 



(58) 



B. High temperatures, n <^T. 

In this case, the low field asymptotics of MR is realized 
at <C T and is given by Eq. (|58p . where one can put 
M = 

,5.,, 



Hence, near the Dirac point resistivity correction scales 
as \/H at if — >■ 0, and increases (for fixed H) with de- 
creasing the temperature. 

Let us now consider larger fields, ^ T. In this case, 
one can use low-energy asymptotics, 

...(e).^L^,..,(e).^^. (60) 



Thermal averaging yields 

(Jxx = 



de 

Snh Jo elT 

9CiC(3) 6^7 T3 



cosh^ (^) cosh^ (^) 



27r 



h el' 



(61) 



a 



xy 



4 In 2 6^7 /xT 

TT h el' 



cosh^ (^) cosh^ 



(62) 



where C, is Riemann zeta function, ^(3) ~ 1.2. From 
Eqs. dSI]) and (|62]) we find 



^Qxx _ 4 £^ 

Pxx(o) ^ 9CiC(3) r3 



sin 2 



9CiC(3) 



2 99 

P- e l 

rp4 



del/T^, for T < < T^/fi; 
C5£*T//i2, for rV/i<e,. 



(63) 



Here C4 = 4/9CiC(3) w 0.54, C5 = 9CiC(3)/16(ln 2)^ « 
0.96. 

The results of calculations are summarized in the 
Fig. 4. 



IV. HALL COEFFICIENT 



Using the equations derived above one can easily cal 
ilate the transverse n 
and the Hall coefficient 



culate the transverse resistivity Qxy — <yxy/{<^tx + ^'xy) 



R 



Qxy 

H ' 



(64) 



Below we discuss the dependence oi R on H and on the 
electron concentration at zero field, 

poc pO 

n= np{e)i'o{e)de — / t'o(e)(ie. (65) 



The second term in the r.h.s. of Eq. (|65|) is the concentra- 
tion of background electrons which compensate the pos- 
itive charge of the donors in the neutrality point. One 
can easily check that n = for fi = 0. 
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Figure 4: Resistivity of graphene for the cases T jj, and 
T ^ /J, as a function of e* oc y/H. 




Figure 5: Dependence of the Hall coefficient, g^y/H, on 
£. oc VH at r > /i. 



and calculate the integral in the limits £* and T. Do- 
ing so we find with the logarithmic precision: cr^y ~ 
/■Kh)[^iel/T^)\n{T /e^). The longitudinal conductiv- 
ity is given by Eq. (j52p where one can neglect small 
correction proportional to thus writing Uxx ~ 

2e^7/7rft. Using these equations we find 



T 

R = Ro In I — 



(70) 



where 



For T ^ simple calculation yields that the Hall 
coefBcient up to small corrections is given by conventional 
expression 



R 



in the whole interval of . In this case, 



(66) 



(67) 



Consider now vicinity of the Dirac point, fi ^ T. In 
this case, 



while the transverse conductivity is given by 



(68) 



de 



AT 



1 



^io ''"-^'^C0sh3(^) 



cosh'(^) cosh2(^) 
sinh(^) 



(69) 



For <g; r, the main contribution to <Txy comes 
from the energy interval e, < e < T, where integral 
in Eq. ()69p is logarithmically divergent. Therefore, one 
may use large-e asymptotic, axy ~ (2e^7/7r?i)(2eJ/e^), 



i?o = 



ecNT^ 



n 

7^ 



(71) 



where Ut ~ NT"^ jv^'f? is the electron concentration for 
H ^ T. Hence, the Hall coefficient logarithmically in- 
creases with decreasing the magnetic field. Above we 
noticed that our calculations are valid up to the ex- 
ponentially small fields where e, ~ Ae"'^'''/^. There- 
fore, the maximal value of the ln(r/£*) is limited by 
7r7/2-ln(A/T). 

In the opposite case, ^ T, the conductivity tensor 
is given by Eqs. (|FT|) and which yield for the Hall 
coefficient the following expression 



R = R, 



1 + C7fi^el/T 



(72) 



where Cg = 641n2/[9CiC(3)]2 « 0.82 and C7 = 
[81n2/9CiC(3)]2 « 0.57. We see that the Hall coefficient 
linearly increases with magnetic field, i? (x oc , in the 
interval T ^ ^ ^'^id saturates when becomes 

larger than From Eqs. ([7T|) and ((72)) we conclude 

that R is non-monotonic function of the magnetic field 
and has a minimum (for positive ^) at magnetic fields cor- 
responding to ~ T. The minimal value of R is given 
by Rmin ~ ^0- The dependence of R on is plotted 
schematically in Fig. 5. 

Measurements of the Hall coefficient are usually used 
for extracting the density of the carriers with the help of 
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conventional expression ((66|) . From Eqs. (|68)) . (|70p . and 
([7^ we see that such a procedure fails near the Dirac 
point. Let us, therefore, discuss the dependence of R 
on n in details (see also Ref. HB) ■ The analysis of above 
equations shows that i? is a non-monotonic function of n 
both at low and high fields: it turns to zero for n = 0, 
has a maximum at certain n = rim and decays as 1/n 
according to Eq. for n — )• oo (see Fig. 6). 



^max 




Figure 6: Dependence of the Hall coefficient on electron 
concentration 



As follows from Eqs. ((681), dZQ]), and at low fields 
or, equivalently, high temperatures (e, <C T), the Hall 
coefficient linearly increases with n at n ^ tlt (/i ^ T), 



R 



ecn% 



ln(r/£,), 



reaches the maximum value. 



R„ 



1 



ecTiT 



(73) 



(74) 



at n„j TiT {fi ^ T)^ and decays inversely proportional 
to the concentration at n — ?> oo.— 

At high field or low temperatures (e, ^ T), the depen- 
dence of i? on ri also has a maximum and is linear at small 
concentration. However, as seen from Eq. ()72|) there are 
some differences compared to low-field case. First of all, 
the coefficient in the linear dependence at small n turns 
out to be different 



R 



ecriT, 



(75) 



Secondly, the maximum value of R is reached at much 
smaller concentration 7i,„ ^ {T/e^)nT corresponding to 
fi ~ T^/e* ^ T. Equation (|75p holds below this concen- 
tration while at larger n the Hall coefficient is given by 
conventional expression (j66|) . The maximum Hall coeffi- 
cient reads 



Rn 



1 



ecriT V T 



(76) 



CHARGED IMPURITIES 



In the previous sections we discussed scattering by the 
short-range potential. One can see from above deriva- 



tions, that the most interesting result, the square-root 
dependence of MR on H at low H is a direct consequence 
of the energy dependence of scattering time specific for 
such a scattering: l/rtr oc \e\. 

Next we discuss scattering by the charged impurities, 
which are often considered to give a dominant contri- 
bution to the resistivity. If we neglect screening of such 
impurities, we get l/rtr oc |e|~^, which implies that uJcTtr 
does not depend on energy and, consequently, the mech- 
anism discussed above does not work. We will see, how- 
ever, that the screening dramatically changes the situa- 
tion and, remarkably, the low-field MR is proportional to 
^/H, in a full analogy with short range scattering, though 
the temperature dependence of MR is different. 

The matrix element of the scattering on a single 
charged impurity is given by 



1 + {2T:e'^ /qx)Nn{q) ' 



(77) 



where n(g) is the static polarization operator and x is 
the dielectric constant. Neglecting screening and using 
golden rule we find that transport scattering rate is in- 
versely proportional to the energy. Let us now take 
screening into account. For simplicity, we will consider 
the low field asymptotic only focusing on the MR at the 
Dirac point (/i ^ T). 

The general expression for polarization operator at the 
Dirac point was derived in Ref. [2^ For our purposes it 
is sufficient to know the asymptotical expression for Il{q) 
at g < T/hv : 

^{q)~^-^, torq^T/hv. (78) 

Indeed, as we discussed in the previous sections the main 
contribution to resistivity at weak fields comes from the 
low energies, e <^ T. Since transferred momentum hq is 
on the order of e/w, we conclude that relevant momenta 
are smaller than T/hv. Assuming also that /hvx ^ 1 
we find that 



1 _ ^h^v'' 

iT iVTln2' 



(79) 



We sec that at low electron energies scattering matrix 
element does not depend on the energy just as in the 
case of the short range potential. Hence, in order to find 
low-field MR one should make the following replacement 
Uo ^ Mg = tt1t?v'^ / NT h\2 and, consequently. 



7^7 



2jV^(ln2)^r^ 



(80) 



Since 7 becomes temperature-dependent we conclude 
that energy e* now also depends on T: 



Low- field MR is still given by Eq. ([59)) , where one should 
substitute expression (j81|) for . It is worth emphasizing 
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that charged impurities yield temperature independent 
MR in contrasts to inverse temperature dependence of 
MR in the case of short range potential. The latter fact 
would allow one to distinguish these two scattering mech- 
anisms in experiments. 

VI. CONCLUSIONS 

To conclude we studied magnetotransport in graphene 
with the short-range disorder within SCBA. We found 
that magnetoresistance depends on three relevant param- 
eters having dimensionality of the energy, fj,,T and 
the field dependence being fully absorbed by which is 
proportional to ^/H. 

One of the main predictions of our model is the square- 
root field dependence of magnetoresistance in the limit of 
low H, both at very low and at very high temperatures. 
Such a dependence persists down to exponentially small 
fields, corresponding to ^ Ae^'^'*'/^. 

We separately analyzed the cases of low and high tem- 
peratures and identified four different transport regimes 
for /X » T and three regimes for <C T (sec Fig. 3). 
All these regimes can be realized provided that temper- 
ature lays within a certain interval. Let us now find the 
corresponding criteria. 

For ^ ^ r and not too small e*, MR is determined 
by energies close to the Fermi surface, £ sa /i. Above 
we assumed that Shubnikov de Haas oscillations are sup- 
pressed by energy averaging within the temperature win- 
dow, which implies that T 3> hucifJ')- Using Eqs. (jS]) and 
pop this inequality can be rewritten as <^ ^^T^. 
While identifying regimes shown at the upper picture in 
the Fig. 3 we implicitly assumed that \J^T\i is larger than 
/X, or, equivalently, T > fj,/^. When T becomes smaller 
than T > fi/^ the high-field square-root asymptotic of 
MR is not realized because of arising of Shubnikov de 
Haas oscillations at ~ VjTj^ ^ However, the low- 



field square-root asymptotic is determined by energies 
e ~ £* and remains valid even at low temperature, be- 
cause Eq. (fT2|) is always satisfied for T <C ^ and e» < £2. 

In the opposite limiting case, T 3> /x, the main con- 
tribution to MR comes from £ w £,. Rewriting Eq. (fT2|) 
as £* <^ 7T, we find that three regimes shown in Fig. 3 
are realized provided that 7T 3> or, equivalently, 

T <^ 7/i. At higher temperatures, T ^ 7/j,, MR is given 
by C3£*/2T for £* < T and C^e^/Tf in the interval 
T ^ £» ^ 7T. The Shubnikov dc Haas oscillations ap- 
pear at £* 3> 7T. 

We also predict a non-trivial behavior of the Hall co- 
efficient on H in the vicinity of the Dirac point. With 
increasing the field, R decreases, reaches a minimum and 
then starts to grow again. Further, we analyzed depen- 
dence of R on electron concentration and found that this 
dependence is non-monotonic both for law and strong 
fields. 

Importantly, the main prediction of our theory, the 
square root dependence of MR in the limit of weak fields, 
is also valid for the case of charged impurity potential. 
It is worth finally noting that our theory is applicable 
for any other 2D electron system with the linear energy 
spectrum such as topological insulators 1^1^ 

Note added: After this work was completed the ex- 
perimental evidence of square-root MR in monolayer 
graphene was reportedi^ 
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/(e) 



ecn{e, H) _ axx{e)z{£) 

H LOc{e)Ttr{e)' 



This is the classical relation between (Jxx{£) and (Jxyi^e) 
where the Ttr is replaced with Ttr/z. Subtracting cr^j,(e) 
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III \ 



z{e')\e'\de' 



\e\ez{e) 
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